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Abstract 



We study delay-independent stability in nonlinear models with a distributed delay 
which have a positive equilibrium. Such models frequently occur in population dynam- 
ics and other applications. In particular, we construct a relevant difference equation 
such that its stability implies stability of the equation with a distributed delay and a 
finite memory. This result is, generally speaking, incorrect for systems with infinite 
memory. If the relevant difference equation is unstable, we describe the general delay- 
independent attracting set and also demonstrate that the equation with a distributed 
delay is stable for small enough delays. 
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X ■ 1 Introduction 



In models of population dynamics which are described by an autonomous differential equation 

^ = f(N)-g(N), (1) 

where f(N) and g(N) are reproduction and mortality rates, respectively, f(N) > 0, g(N) > 
for N > and f(N) > g(N) for < N < K, f(N) < g(N) for iV > K (K is the 
carrying capacity of the environment), the positive equilibrium K is stable: all positive 
solutions converge to K and are monotone. It was argued that the observed data usually 
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oscillates about the carrying capacity; in order to model this phenomenon, it was suggested 
to introduce delay in the production term 

dN 

— = f{N{t-T))-g{N)- (2) 

the latter equation can have oscillatory solutions, and the delay incorporated in the right 
hand side can be interpreted as maturation, production or digestion effects. It is usually 
assumed that the mortality rate was proportional to the present population level 

dN 

a — = -^N{t) + f{N{t-r)) 1 M >0. (3) 

The global behavior of solutions of (j3J) has been extensively studied in literature, in particular 
in the cases of negative and positive feedback (see, for example, [D,[2] and references therein), 
the chaotic behavior is impossible in the case of the monotone feedback [3] . However in the 
case when f(x) is a unimodal function, i.e., increases for x < K and decreases for x > K, 
there may be delay induced instability and complex dynamics [HE]. For a detailed overview 
of the literature on the dynamics of see the recent papers O [7] . It is demonstrated in [6j [7J 
that if / is a unimodal function and positive equilibrium K of the equation x n+ \ = f(x n ) 
is globally asymptotically stable, then all solutions of (j3J) tend to K. In particular, if / 
has a negative Schwarzian derivative, then local stability of the equilibrium of the difference 
equation implies its global attractivity [5]. To the best of our knowledge, the first delay- 
independent stability conditions were obtained in [9]. In the present paper we will try to 
answer the general question: what are intrinsic properties of the reproduction function / 
which allow us to conclude that any solution of the equation with a finite memory converges 
to the equilibrium? Here we consider both general delays (including integral terms) and 
continuous functions / which may have multiple extrema, tend to infinity at infinity etc. 

As special cases, (EJ) includes the Nicholson's blowflies equation [101 [11] and the Mackey- 
Glass equation [U [13] . The Nicholson's blowflies equation 

x(t) = -5x(t) + px{t - T)e- ax(t - T) (4) 

was used in [11] to describe the periodic oscillation in Nicholson's classic experiments [10] 
with the Australian sheep blowfly, Lucila cuprina. Equation (j4j) with a distributed delay 
was studied in [T2], where comprehensive results were obtained for the case 5 < p < 5e. 
The Mackey-Glass equation [U [T3] 

= ax(t-r) _ 
w 1 + afr(t - r) w w 

models white blood cells production. Local and global stability of the positive equilibrium 
for equation ([5]) with variable delays was studied in [HI [TJ2, [161 E0 EEl EH]; to the best of 
our knowledge, there are no publications on (jSJ) with a distributed delay. 

To incorporate random environment influence, some authors included noise in (j3J) and 
studied attractivity conditions, see, for example, [20]. However, in applied problems not only 
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the derivative but also the delay value can be perturbed. We assume that the production 
delay is not a constant r but some distributed value which leads to the equation 



x(t) 



f(x(s))d s R(t, s)-x{t) 



(6) 



t-b 



where / d s R(t, s) is the probability that at time t the maturation delay in the production 



t—a 



function is between b and a, where < b < a. We will assume that very large delays are 
improbable, substituting — oo in the lower bound with h(t) < t which tends to infinity as 
t — > oo. In the present paper we consider a rather general form of /, which includes unimodal 
functions, as well as functions with several extrema. The only requirement is that f(x) has 
the only positive fixed point. The main result claims that if this fixed point is a global 
attractor for all positive solutions of the difference equation 



%n+l — f{ x n) 



(7) 



then all solutions of ([6]) with positive initial conditions tend to this fixed point as well. To 
some extent this establishes a link between stable differential equation (TjQ) and difference 
equation ([7]) which can undergo a series of bifurcations and even transition to chaos. If ([7]) 
is globally stable, so is If the unique positive equilibrium of ([7]) is unstable, §6§ can be 
stable or not, depending on the delay. 

The paper is organized as follows. In Section 2 we prove that all solutions with posi- 
tive initial conditions are positive and bounded and establish some estimates for the lower 
and the upper bounds. Section 3 presents sufficient conditions under which all positive 
solutions converge to the positive equilibrium. In Section 4 delay-dependent stability is in- 
vestigated. In particular, it is demonstrated that equations are globally attractive for delays 
small enough; if (JTj) is unstable, then we can find such delays that the positive equilibrium of 
([6]) is not a global attractor. In Section 5 these results are applied to equations of population 
dynamics with a unimodal reproduction function and a distributed delay, in particular, to 
the Nicholson's blowflies and Mackey- Glass equations; some open problems are presented. 



2 Boundedness and Estimates of Solutions 

We consider the equation with a distributed delay 



x{t) = r(t) 
and the initial condition 
As special cases, OH]) includes 



f(x(s))d 8 R(t,s)-x(t) 



Uhit) 



t > 0, 



x(t) = <p(t), t < 0. 



(9) 



3 



1. The integrodifferential equation 



x(t) = r(t) 



K(t, s)f(x(s))ds-x(t) 



Uh(t) 



(10) 



corresponding to the absolutely continuous R(t, •) for any t. Here 

r* d 

/ Kit, s) ds = 1 for any t, Kit, s) = -^-R(t, s) > 
Jh(t) ds 

is defined almost everywhere. 
2. The equation with several concentrated delays 

m 

J2a k (t)f(x(h k (t)))-x(t) 



x(t) = r(t) 



k=l 



with a k {t) > 0, k = 1, • • • ,m, where Qfe(^) = 1 for any t. This corresponds to 

fc=i 

s) = cik{t)x{h k (t),oa) (s) , where Xi{t) is the characteristic function of interval /. 
fc=i 

Definition. An absolutely continuous in [0, oo) function x : R — > R is called a solution of 
the problem (ED,® if it satisfies equation (jSJ) for almost all t G [0, oo) and conditions ([9]) for 
t < 0. 

The integral in the right hand side of (jHJ) should exist almost everywhere. In particular, for 
(TTTJT) with a locally integrable kernel, (p can be any Lebesgue measurable essentially bounded 
function. For (llip 9? should be a Borel measurable bounded function. For any distribution R 
the integral exists if if is bounded and continuous (here we assume / is continuous). Besides, 
as is commonly set in population dynamics models, (p(t) is nonnegative and the value at the 
initial point is positive. 

Consider (JED,© under the following assumptions. 

(al) / : [0,oo) — > [0,oo) is a continuous function satisfying Lipschitz condition \f(x) — 
f(y)\ < L\x - y\, x,y > 0, /(0) = 0, f(x) > x for < x < K and < f(x) < x for 
x > K; 

(a2) h : [0, 00) — > R, is a Lebesgue measurable function, h(t) < t, lim hit) = 00; 

t— »oo 

(a3) r(t) is a Lebesgue measurable essentially bounded on [0, 00) function, r(t) > for any 



t > 0, 



ris) ds = 00; 



(a4) R(t, •) is a left continuous nondecreasing function for any t, R(-, s) is locally integrable 
for any s, R(t,s) = 0, s < h(t), R(t,t + ) = 1. Here u(t + ) is the right side limit of 
function u at point t. 
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(a5) if : (—00, 0] — > R is a continuous bounded function, ip{t) > 0, y?(0) > 0. 

First, let us justify that the solution of (jHJ),© exists and is unique. 

Denote by L 2 ([t ,ti]) the space of Lebesgue measurable functions x(t) such that 




by C([io,ii]) the space of continuous in [to,t\] functions with the sup-norm. 

We will use the following result from the book of Corduneanu [21] (Theorem 4.5, p. 95). 
We recall that operator N is causal (or Volterra) if for any two functions x and y and each 
t the fact that x(s) = y(s), s <t, implies (Nx)(s) = (Ny)(s), s <t. 

Lemma 1 [21] Consider the equation 

y(t) = (jCy)(t) + (Afy)(t), te Mi], (12) 

where £ is a linear bounded causal operator, N : C([£o,£i]) —* L 2 ([to,£i]) is a nonlinear 
causal operator which satisfies 

\\Afx - Afy\\tf(it , tl ]) < X\\x - Z/||c([t ,tiD ( 13 ) 

for A sufficiently small. Then there exists a unique absolutely continuous solution of (Tjj|) in 
[to,ti], with the initial function being equal to zero fort < t Q . 

Theorem 1 Suppose (al)-(a5) hold. Then there exists a unique solution of 

Proof. To reduce ([8]) to the equation with the zero initial function, for any ( > we can 
present the integral as a sum of two integrals 

x(t) = -r(t)x(t) + r(t) [ f(x(s)) d s R(t, s) + r(t) [ f(<p(s)) d.R(t,s), (14) 

J t J t 

where 

x{t) = 0, t < t , (p(t) = 0, t > h. 

Here to > is arbitrary, so we begin with to = and proceed to a neighboring t\ to prove 
the existence of a local solution. Then in (fl~2l) 

(£x)(t) = -r(t)x(t), (Afx)(t) =r(t) [ f(x(s))d s R(t,s)+F(t), 

Jto 

where 

F(t) = r{t) [ f{<p(s)) d s R(t, s), x{t) = 0, t < t , <p{t) = 0, t > t , 

Jt 
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and for any A > there is ti, such that 



\f(x(s))-f(y(s))\d s R(t,s) 



to 



< Less sup \r(t)\ max \x(s) — y(s 



e[to,*i] 



L 2 ([*o,*i]) 



d s R(t, s) 



t 



L 2 ([t ,ti]) 



< Lesssup \r(t)\ \\x(s) - 2/(s)||c([to,ti])l*o - *i| < A||ar - y||c([*d,*il) 



for | — ^1 1 < A/(Lesssup |r(t)|), where L was defined in (al), here A can be chosen small 



t>o 



enough. By Lemma [T] this implies existence and uniqueness of a local solution for (jSJ). This 
solution is either global or there exists £2 such that either 



or 



liminf x(t) = —00 
t-»t 2 



limsup x(t) =00 . 
t->t 2 



(15) 



(16) 



The initial value is positive, so as far as x(t) > 0, the solution is not less than the solution 
of the initial value problem x + r(t)x = 0, x(0) = Xq > which is positive and the former 
case ( TTB1) is impossible. In addition, x(t) < for any 



x(t) > K, x(t) > max < max x(s), sup ip(s) > , 

[0<s<t s < J 

which contradicts (TT6"|) . Thus there exists a unique global solution, which completes the 
proof. □ 

Theorem 2 Suppose (al)-(a5) hold. Then the solution of (E]) ; (EJ) is positive for t > 0. 
Proof. After the substitution 



y(t)=x(t)ex P \ / r(C) 



(17) 



equation (jSJ) becomes 

y(£) = r(t) exp < / r(s) (is 



/ h/(s) exp 



r(C)dQ t>0. (18) 



h(t) 

Thus ?/(0) > and y(t) > as far as y(s) > 0, s < t, consequently, y(t) > for any t > 0. 
Since the signs of and x(t) coincide, then x(t) > for any t > 0. □ 

Definition. The solution x(£) of (151), (151) is permanent if there exist A and 5, B > A > 0, 
such that 

A < x(t) <B, t>0. 

In the following we prove permanence of all solutions of (|HJ) with positive initial conditions; 
moreover, we establish bounds for solutions. 
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Theorem 3 Suppose (al)-(a5) hold. Then a solution of is permanent. 



Proof. By Theorem [2] the solution is positive for t > 0. By (a2) there exists to > such that 
h(t) > 0, t > to. Since the solution is a continuous positive function, then we can define 



x min = min x(t) > 0, x 
te[o,t ] 



max — max ./'(/). 

te[o,t ] 



(19) 



Without loss of generality we assume x min < K, x max > K; otherwise, we can choose 
min{x min , XK}, max{x max , K/X}, where < A < 1, as x min and x max , respectively By (al) 
the following values are positive 



Define 



M = max f{x), m = min f(x) 

3Jl£: [x mi n ; 2^- max ] [^mini^max] 



B = max < M, x max , max f(x) > , A = min < m, x m i n , min f(x 



xe[o,K] ' 



xe[K,B] ' 



(20) 



(21) 



Since f(x) > x, < x < A and f(x) < x, x > B, then there exists 5 > such that f(x) > A 
for A — 5 < x < B and /(x) < i? for A < x < B + 5. Let us demonstrate 



x(t) E [A,B], t> 0. 



(22) 



By the definition of A, B we have x(to) G [A, B\. Suppose the contrary: x(t) > B or x(t) < A 
for some t > to- 

First, let a:(t) > B for some t > to- Then x(t) = B + e for some e < 5. Denote 

S 1 = {t>t \x(t)> B + e}, t* = infSi. 
Since x(to) < B then the set 

S 2 = {t\t <t<t*, x(t) < B} 

is nonempty, denote t* = sup S2. Then x(t*) = 5, x(t*) = B + e; we also have B < x(t) < 
B + 8 and /(x) < 5 in the interval [t*,t*], thus the derivative is nonpositive 



i(t) = r(t) 



/(x(s)) d s R(t,s)-x(t) 



Uh(t) 



< r(t) 



B d a R(t, s)-B 



Uh(t) 



0. 



which contradicts the assumption x(t*) = B + e > B = x(t*). 

Similarly, let us assume that x(t) = A — e for some e > 0, e < 5 and some t > to- After 
introducing 

S 1 = {t > t \x(t) < A - e}, t* = inf S h S 2 = {t\t <t<t*, x{t) > A} } U = sup S 2 , 
we have x(Q = A, x(t*) = A - 5, A - 5 < x(t) < A and f{x(t)) > A for t E [U,t% hence 



x(t) = r(t) 



f(x(s)) d s R{t,s)-x(t) 



h(t) 



> r(t) 



/ A d s R{t, s)-A 

Jh(t) 



0. 



This contradicts the assumption x(t*) < x(t*). Consequently fl22l) is valid for t > t and 
thus for any t > 0, the bounds are positive, so the solution is permanent, which completes 
the proof. □ 

Example 1. The statement of Theorem [3] is not valid if we omit the condition lim h(t) = oo. 
Consider the equation 

x(t) = 5x(h(t))e- x(h(t)) - x(t), t > 0, h{t) = -1, x(t) =t + l, t G [— 1, 0], (23) 

which is equivalent to the initial value problem 

x(t)+x(t) = 0, x(0) = 1, (24) 

its solution x(t) = e~* tends to zero as t — > oo and so is not permanent. 



3 Absolute Global stability for Stable Difference Equa- 
tions 

One of the main steps in establishing global stability property is the proof of the fact that 
all nonoscillatory about the equilibrium solutions tend to this equilibrium (see, for exam- 
ple, [22\). For ordinary differential equations all solutions are nonoscillatory, for retarded 
equations it depends on the delay. Below we demonstrate that convergence of nonoscillatory 
solutions to the equilibrium is quite a common property which is valid for any reproduction 
function with a unique positive equilibrium. It can be interpreted as: "if nonoscillatory, so- 
lutions of delay equations behave asymptotically similar to ordinary differential equations" . 



Definition. A solution x(t) of (IE]),© is nonoscillatory about K if there exists r > such 
that either x(t) > K or x(t) < K for all t > r. Otherwise, x(t) oscillates about K. 

Theorem 4 Suppose (al)-(a5) hold. Any nonoscillatory about K solution of (Et^iE)) con- 
verges to K . 

Proof. First, let x(t) < K, t > r. Without loss of generality we can assume r = 0. By (a2) 

there exists to > such that hit) > for t > to- Denote A as in (12TI) . By Theorem [3] we 

obtain that x(t) > A for any t > 0. Since / is continuous and f(x) > x for x < K, then 

m = inf f(x) > A. There may be two cases: m = K and m < K. In the former case, 
xe[A,K] 

since x(t) < K, we have 



x(t) > r(t) 



xE[A,K] 



> 0, t>t, 



li- 



as far as x(t) > A, thus the solution of the delay differential equation is not less than the 
solution of x(t) = r(t)[K — x(t)}, < x(to) < K, which is increasing and by (a3) (the integral 
of r(t) diverges) tends to K. 
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A lfl m,) I, m, K=M, M, M B x 

Figure 1: For an arbitrary reproduction function with one positive equilibrium K we con- 
struct a series of such points that eventually a solution is in [rrij, K], if it does not exceed K 
and is in [K, Mj] if a solution is not less than K. 



Consider the latter case m < K. By the definition of mo and f(x) > x for x G [A, K] we 
have lo = sup{x < K\f(x) < m } < m . Taking any a, lo < a < m and assuming x(t) < a 
for any t, we obtain 



x(t) = r(t) 



f(x(s)) d s R{t,s)-x{t) 



Uh(t) 



> r(t)(m - a) > 0, 



which leads to a contradiction x(t) — > oo as t — » oo since f °° r(s) ds diverges. Thus, 
> a for some t*; moreover, since x(t) > as a;(t) < m then x(t) > a for any t > t*. 
Let /t(t) >t*,t>t* for some t*. By the definition of lo and a we have m = inf /(a;) > mo 

and as far as x(t) < mo and h(t) > t* the following inequality holds 



x(t) = r(t) 



f(x(s)) d s R{t,s)-x{t) 



Uh(t) 



> r(t)(m — m ) > 0. 



Assuming x(t) < mo for any t we again obtain a contradiction. Thus, there exist fii and 
t\ > t such that x(fii) > m and h(t) > \ii for t > t\. Then > m for any t > jii and 
x(h(t)) > m Q , t > ti. 

Further, let mi = inf f(x) < K, here mo < mi since /(x) > x for < x < K. 

x&[mo,K] 

Similarly, we obtain x(t) > mi whenever t > £2, for some t<i > t\. 

We continue this process. It can be finite (for example, in Fig. [T]we have m2 = K, where 
the process stops and we deduce x(t) — > K as t — > 00) or infinite (see the branch x{t) > K 
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of Fig. 1). In the infinite case we have an increasing sequence {rrij}, rrij + i = min f(x) 

x£[rrij,K] 

which does not exceed K, so this sequence has a limit d. Since f(x) is continuous then 
d = min fix). If d < K then f(x)—d should attain its minimum in [d, K\ but fix) > x > d, 

x£[d,K] 

so this minimum is positive and the equality d = min fix) leads to a contradiction. 

xe[d,K] 

Further, let x(t) > K. Similarly, we define B as in (121!) and Mq = max f{x) < B. There 

x€[K,B] 

may be two cases: Mq = K and Mq > K. In the former case we obtain x(t) — > K as t — > oo. 
Consider the latter case. By Theorem[3]we have x(t) < B for any x > 0. By the definition of 
M and f(x) < x, x G [K, B] we have so = inf{x > K\f(x) > M } > Mq. Similar to the case 
x(t) < K we demonstrate that there exists V\ such that xiv\) < Mq and x(t) < Mo for any 
t > V\. Let ti be such that hit) > V\ for t >t\. Further, we define Mi = max fix) < Mq. 

xe[K,Mo] 

We continue this process, it can be finite or infinite (Fig. [1] illustrates an infinite process for 
x(t) > K). Similar to the case x(t) < K we obtain x(t) — >• K as t — > oo, which completes 
the proof. □ 

Example 2. Let us note that in the case of infinite delays nonoscillatory solutions do not 
necessarily tend to the positive equilibrium. For example, the solution of the equation 

x (t) + x {t) = | x(0) e~ x{0 \ x(0) = 2, 



which is xit) = — h ( 2 e , tends to lie while the positive equilibrium is 1 — In 2, the 

e \ ej 

monotone solution is nonoscillatory. 

Thus for any reproduction function with a unique fixed point fix) = x nonoscillatory 
solutions tend to the equilibrium; this is not, generally, true for oscillatory solutions. 

Example 3. Consider the Nicholson's blowflies equation 

x(t) = -6x(t) +px(t - r)e- ax{t - T) . (25) 

Denote 



1 



5^aN*(aN* -2) 



. JaN*(aN* - 2. . 

arcsm — + 2nk 

aN* — 1 



fc = 0,l,2,---, (26) 



where iV* = l/a\n(p/8) is a positive equilibrium. If p > 5e 2 then the positive equilibrium is 
locally asymptotically stable for r € [0, r ) and is unstable (locally, thus it cannot be globally 
attractive) for r > r , and ( |25i) undergoes a Hopf bifurcation at iV* when r = [23] for 
k = 0,1,2, 

Now we prove that absolute (delay-independent) convergence holds in some special cases. 
Lemma 2 Suppose (al)-(a5) and at least one of the following conditions holds: 
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1. f(x) < K for any < x < K . 



2. Denote by x max ([a, K]) , < a < K, the greatest point in [a,K] where max f(x) is 

xe[a,K] 

attained and assume 

f (x ma J\a, K})) > K for some a G (0, K) =>- min fix) > a. (27) 

X£[K,f(x max ([a,K])]y 

Then any solution of ^),^) converges to K. 

Proof. 1. By Theorem [3] any solution is permanent: A < x(t) < B, where A < K < B, 
t > 0, and for some t > we have h(t) > 0, t > t . Denote M = max f(x) and 

x£[K,B] 

inn = min f fx). 

xe[A,K] V ' 

If M < K then the derivative is negative for any x(t) > K and the solution either 
eventually does not exceed K or is decreasing for any t > to. In the latter case the solution 
tends to the equilibrium; if it has a different limit, we obtain that the derivative is less than a 
negative number, which is a contradiction. In the former case, if there exists t > such that 
x(t) < K (otherwise, we have a nonoscillatory case x(t) > K where we have already proved 
convergence), then x(t) < K for any t (assuming x(t*) > K we obtain that the derivative of 
x(t) is negative almost everywhere, while the function changes from K to x(t*) > K), which 
is again a nonoscillatory case, by Theorem H] solution x(t) converges to K. 

Thus, we can consider Mq > K only. Then we introduce 

Mi = max f(x) < Mo, mi = min > tuq, Si = min f(x) 

xe[K,M ] x£[m ,K] x£[K,M ] 

and similarly define Mj, mj and Sj, j = 2, 3, • • • . There exists T\ such that min{mi, s±} < 
x(t) < Mi for t > T\ and t\ such that h(t) > r%, t > t\. Similar to the proof of Theorem H] 
we obtain that there exists a sequence t± < t% < • • • < Tj < • • • such that 

mm{rrij, Sj} < x(t) < Mj for t>Tj. 

We assume that all Mj > K, otherwise we have an eventually monotone case. Since for any 
continuous / all three sequences are monotone (nonincreasing Mj and nondecreasing rrij, Sj) 
and tend to K (see the end of the proof of Theorem H]), then lim x(t) = K. 

t^oo 

2. Now suppose that fl27j) holds for any a > 0. We begin with A < x(t) < B in Theorem 
El t > 0, where A, B are defined in (J2D, h(t) > 0, t > t . Denote 

M = max fix), m = min f(x). (28) 

xe[A,K] xe[K,M ] 

The case M = K was considered in Part 1, thus we can restrict ourselves to the case 
Mo > K, maximum is attained at 

xq = x max _([A,K]), 
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here Xq is the greatest point where the maximum is attained. Let us demonstrate that there 
exists fjL such that x < x(t) < M , t > fi . 

If x(t) > Sq = max{M , max fix)} = max fix) > K for all t then we have a nonoscil- 

x£[K,B] xe[A,B] 

latory case and convergence to K, which is a contradiction. Thus, x(t*) < Sq for some 
t*; assuming there is t* > t* such that ^(t*) > So we obtain that the value of the func- 
tion at the end of the segment is higher than at the beginning point, while the derivative 
is nonpositive. So there exists r > to such that x(t) < So for t > r . If Sq ^ M then 

Si = max f(x) < So. Similarly, we find t± such that x(t) < Si for t large enough. Since 

xe[K,s ] 

the sequence S n +\ — max f(x) is nonincreasing and tends to K then there exists Ti such 

xE[K,S n ] 

that x(t) < M , t > T\. Further, we will consider t > ti, where hit) > t\ whenever t > ti, 
only. 

If x(t) < so = minjmo, min f(x)} = min f(x) for all t, where mo was defined in (128|) . 

xG[A,K] xe[A,M ] 

then we have a nonoscillatory case and convergence to K, which is a contradiction. As above, 
we prove that x(t) > s for t > u , here u > t . If s < ttiq, then we construct a sequence 
s n+ i = min f(x) which tends to K. If mo < K then some Sj > mo- There is fio > t 1 such 

x£[s„,K] 

that x(t) > m , t > fio. Consequently, we have found fio such that m < x(t) < M , t > fio. 
Let h{t) > fio, t > t 2 . 

Now, if mo = K then for t > t 2 the solution increases as far as x(t) < K. Thus, either 
x(t) < K for any t > t 2 and this monotone solution converges to K, or, if x(t*) > K for 
some t*, x(t) > K, t > t*, and again we have a nonoscillatory solution which converges to 
K. 

Denote 

Mi = max f{x), mi = min f(x). 

xe[mo,Mo] x£[mo,Mi] 

Let us assume Mi > K, mi < K and demonstrate that there is fii > t 2 such that mi < 
x(t) < Mi, x(t) > fii, where Mi = f (x max ([m , K])). In fact, for any t > t 2 the solution 
is nonincreasing as far as x(t) > Mi, which gives an upper bound. Considering t where the 
equation refers only to the values where this bound is valid, we obtain that the solution is 
nondecreasing if x(t) < mi, which together with mi < K < Mi confirms the statement. We 
continue the induction process 

M n+ i = max f(x), m n+1 = min f(x), 

x<={m„,M„] x£[m n ,M„ + i] 

where K > m n+ i > m n , K < M n+ i < M n . This process can be infinitely continued if all 
m n < K, M n > K (otherwise, at certain stage we have a "monotone" case which implies 
convergence), and there exists fij, j — 0, 1, 2, • • ■ such that mj < x(t) < Mj, t > fij. Let us 
assume that there are infinite sequences {m n } and {M n }, both are monotone and bounded. 
Then there exist limits 

x = lim mj < K < lim Mj = X. 
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Since by the assumptions of the theorem f(Mj) > rrij, then f(X) > x. If x = f(X) then 
either x = X = f(X) = X and x(t) converges to K or min f(y) < x which 

V£[K,fmax([x,K])] 

contradicts the assumptions of the theorem. 

Denote by x max j, x m inj the sequences where minimum rrij and maximum Mj are attained, 
respectively (we recall that we choose maximal x maxJ - and minimal Xminj if an extremum is 
attained at several points). Then /(x maxJ ) = Mj, f(x min j) = rrij, rrij < x max j < K, K < 
x m i n ,j < m j, m j+i > Xmax,j, < x min j, which implies that x max j and x min j tend to x and 

X, respectively. Let us prove that x = f{X). Assume the contrary: f(X) > x+2e > 0. Since 
lim ic min „• = X and / is continuous then there exist n such that f(x min »•) > x + e, j > n , 

but /(x m inj) = rrij+i, so rrij+x > x + e, which contradicts the equality lim rrij = supm^ = x. 
Thus lim x(t) = K, which completes the proof. □ 
Now we can prove the main result of the present paper. 

Theorem 5 Suppose (al)-(a5) hold and all positive solutions of the difference equation 

x n+1 = f(x n ) (29) 

tend to K. Then any positive solution of (TJj), converges to K. 

Proof. If either conditions of Part 1 of Lemma [2] or (|27|) hold then K attracts all positive 
solutions of Let us assume that (I2~T1) does not hold (and we do not have a monotone 

case as in Part 1 of Lemma [2]), which means that for some a G (0, K) we have 

min f{x) < a, or K <b < /(a), fib) < a, where b e [K,f(a)}. 

x£[KJ{a)} 

We can assume f(b) = a, otherwise, since f(b) < a, f(K) = K > a and / is continuous, 
then there is c G [K, b) such that f(c) = a. Thus, b > K and f 2 (b) = f(f(b)) = f(a) > b. 
Since fix) < x for any x > K then for any x > M := max fix) we have f 2 (x) < x, then 

xE[0,K] 

there is a fixed point of f 2 in the segment [b, M], in addition to a fixed point x = K. Since 
the fixed point K of (1291) cannot be a global attractor unless K is the only fixed point of f 2 
[21] then K is not a global attractor of (|29[) . which completes the proof. □ 



Here we have not considered the case when / has no positive equilibria; however for 
completeness we will consider this case as well. 

Theorem 6 Suppose (a2)-(a5) hold, /(0) = and f{x) < x for any x > 0. Then any 
positive solution of (t2P,(Tj|) converges to zero. 

Proof. Let to be such that h(t) > for t > to- Denote M = sup x{t). Since f{x) < 

te[o,to] 

x then x{t) < Mq for any t > 0. The solution is decreasing as far as x > M\, where 
Mi = max fix) < M . 

x€[0,Mo] 
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f(x)A 




Figure 2: If f(x) < x for any x > then all positive solutions tend to zero. 



Suppose 5*i is the smallest point not exceeding Mo where this maximum is attained (see 
Fig. [2]). Since f(x) < x, then Si > My. Let us choose li such that M\ < l\ < S\\ as 
far as x(t) > li we have x(t) < r(t)(Mi — h), thus eventually any solution is less than li. 
Further, for some d > we have Mi — f(x) > d if < x < l\. Thus the derivative at any 
point x(t) > Mi does not exceed (we recall that x(t) < h) the negative value of —dr(t). 
Consequently, there exists ti such that x(t) < Mi for t > t±. By induction, we denote 



Mi 



k+l 



< M k , k 



1,2, 



sup f(x 

x£\0,M k ] 

x(t) < Mk+i for t > tk+i- Since lim M n = then lim x(t) 
tends to zero. 



and prove that there exists tk+i > tk such that 

0, i.e., any positive solution 

□ 



4 The Case When Stability Is Delay-Dependent 

In this section we consider the case when stability properties of (JS]) depend on the delay. 
First, we prove that for small delays all solutions tend to positive equilibrium K. Second, 
we demonstrate that if fl27j) does not hold then there exists equation (jSJ) with parameters 
satisfying (al)-(a5) such that its attracting set is as close to [m, M] as prescribed, where 

M = max f(x) > K, m = min f(x) (30) 
xe[o,K1 xe[K,M] y ' y ' 

and max f(x) = M as well. 

x£[m,K] 

Let us note that the Lipschitz condition in (al) implies \f(x) — K\ < L\x — K\. 
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Theorem 7 Suppose (al)-(a5) hold and 



limsup / r(s) ds < 

Jh(t) 



t—*oo 



L + l 



(31) 



Then any positive solution of $Bj),$Q) converges to K. 



Proof. Without loss of generality we can assume that (1311) is satisfied for t > 0. We can 

/"* A 

find a positive A < 1 such that / r(s) ds < — for t > 0. By Theorem [3] any solution 

A(t) -tv + 1 

is permanent: 

# - a < x{t) < K + b , t > 0, < a < K, b > 0. 

Let us denote a = max{a , b } and prove that there exists t\ > such that x(t) < K + Xa 
for any t > t\. If solution x(t) is nonoscillatory about K then by Theorem H] it tends to the 
equilibrium and thus for e = Xa > there exists ti such that x(t) < K + e for t >t\. Now 
assume that is oscillating, t is such a point that h{t) > for t > t and = i^, 

a^r*) = f^, where > for t G (i*, r*). A continuous function x(t) attains its maximum 

in [i*, r*], let = max x(s) = M . We claim that x(t*) < K + Xa. Assume the contrary: 

se[t*,T„] 

x(t*) > K + Aa. By (al) we have max f(x) = M 1 < M . Denote 

x£[K,M ] 

M 2 = max{M 1 ,K + Xa}, r = sup {t G [t*, t* ]\x(t) < M 2 } . 

Then x(t) > M 2 > K + Xa, t G [r, £*]. 

Since — a < a;(i) < K + a for any t > then by (al) 



= r(t) 



/(z(s)) d s R{t,s)-x(t) 



h(t) 



< r(t) 



max f(x) — (K — a) 

x£[K-a,K+a],x>0 



r(t) 



max [f{%) — K] + a 

x£[K-a,K+a],x>0 



< r(t)[La + a] = a(L + l)r(t). 



(32) 



We remark that x(t) < whenever the following conditions hold: t G (r, £*), /i(i) > (all 
referred prehistory of the solution is between K and Ml < M ) and x(t) > M 2 > M\. Since 
/(t) = M 2 , /(£*) > M 2 , then there are points in (r, t*) where the derivative is positive, thus 
h{t) < t*. Let i be such a point. As r{t) > and [t*,r] C then 



A 



r(s) ds < r(s) ds < 
u Jh(t) L + l 



consequently, 

x{t) — x(t*) = / x(s) ds < [ a{L + l)r(s) < a{L + 1) ^ = Aa. 
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Hence x(t) < x(t*) + Xa = K + Xa, which contradicts the assumption x(t) > K + Xa. If 
we denote t\ = then x(t) < K + Xa for any t > t\. 

Similar to the previous case, we can prove that there exists T\ such that x(t) > K — Xa 
for t > t\. 

Thus, we accept K — Xa < x(t) < K + Xa as new solution bounds and consider t 2 such 
that h(t) > max{t 1; ri} for t > r 2 . We repeat the induction step and obtain that there exists 
t 2 > r 2 such that 

K - X 2 a < x(t) < K + X 2 a, t>t 2 , 
and a sequence of points ti < t 2 < t 3 < t n < ■ ■ ■ such that 

K - X 3 a < x(t) < K + X 3 a, t > t 3 , ■ ■ ■ , K - X n a < x(t) < K + X n a, t>t n ,---. 

Since < A < 1 then lim x(t) = K, which completes the proof. □ 

t— »oo 

Corollary 1 Suppose (al)-(a2),(a4)-(a5) hold, r > and sup(t — h(t)) < —— — - (we 

t>o r \L + 1) 

recall that L is the Lipschitz constant defined in (al)). Then all solutions of the equation 



x{t) = r 

tend to K. 



f(x(s))d s R(t, s)-x(t) 
Uh(t) 



(33) 



Remark 1 Let us remark that the Lipschitz constant L in l[31\) can be substituted by, gen- 
erally, a smaller value of 

f(x)-K 



sup 

x>0 



x-K 

since in [3^) only \f(x) — K\ is estimated for \x — K\ < a. 



(34) 



Next, let us estimate the attracting set when ([H]) is not absolutely stable. We introduce 
m, M as in (130]) and denote by x max the greatest point in [a, K] where M = max f(x) is 

xG[0,K] 

attained, by x m j n the minimal point where m = min fix) is attained. 

x£[K,M] 

Theorem 8 Suppose f(x) satisfies (al) and 

m = f(x min ) < x max < K. (35) 
Then for any a G (m, x max ) ; and any b G (K,M) such that min f(x) < a there exists a 

xe[K,b] 

problem (G|) ; (rj|) with parameters satisfying (a2)-(a5) such that 

liminf x{t) = a, limsupa;(t) = b. (36) 



t^oo 
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f(x) A 




Figure 3: Assuming a and b as in the figure exist, we can construct a delay equation such 
that [a, b] is its attracting set. 



Proof. Let us fix a G (m, x max ), and b G (K, M) such that d = min fix) < a, see Fig. |3j 

xe[K,b] 

Then there exists x\ such that K < x± < b and mi = f{x{) < a. For the initial function 

(p(t)=a + (t+l)(b-a), £g[-1,0], 

we have <£>(— 1) = a, <£>(0) = 6 and ip changes continuously from a to b. Since a < K < X\ < b, 
then x(sq) = X\ for some s G (—1,0). We consider the equation 

x(t)=r[f(x(h(t)))-x(t)\, t>0, (37) 

where r > 0. Obviously parameters of f|3T|) and the initial function satisfy (a3)-(a5). 

We choose = so, t G [0, ri], where ri = - In I > 0. Then 



r \ a — m 

Ft 



x(t) = (b - m)e + m, t G [0,Ti], 

so x(ri) = a and x(t) changes continuously from b to a in [0, t±\, a < x max < b, thus 

there exists si G (0, Ti) such that x(si) = x max . We assume /i(t) = s 1; t G (t 1 ,t 2 ], where 

1, [M-a\ 
t~2 = Ti H — In — — - > n. The solution is 



r \M-b 



x(t) = (a - M)e- r{t ~ Tl) + M, £ G [ri, r 2 



hence a;(r 2 ) = 6. We continue periodically with fo(t) piecewise constant such that 
x(h(t)) = { Xl > l^^r 2n+l ] n = 0jlj2 ... j 

^ ^max * ^ ( r 2n+l? r 2n+2ji 

where r = 0, h(t) G (r fc _i, r fe ) if t G (r fc , r fc+1 ), x(r 2n ) = b, s(r 2n+1 ) = a, n = 0, 1, 2 • • • and 
a < x(t) < b for any t. Here the delay is bounded and piecewise continuous, it obviously 
satisfies (a2). □ 
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Corollary 2 If (al)-(a5) and (ESP hold then the sharpest attracting interval for (G|),(G|) is 
(m,M), where m and M are defined in / fffOj) . 

Remark 2 Suppose the difference equation x n+ \ = f(x n ) has a stable 2-cycle (a, b) in such 
a way that f(f(x)) < x for a<x<K,x>b and f{f{x)) > x for x < a, K < x < b. Then 
for any e > there exist such distributed delay and initial function that [a + e, b — e] is an 
attracting interval of |2P ; ([2)]. Moreover, if K e (a, b), function f satisfies (al), /(a) > b 
and f{b) < a, then, similar to the proof of Theorem^ we can demonstrate that there exists 
problem $Bj),$Q) such that (a, b) is its attracting interval. 

Remark 3 The claims of Corollary [H and Remark [H to some extent complement Theorem 
6 in J?|/ where the sharpest invariant and attracting interval was found in the case when 
absolute stability fails but there exists a unique globally attractive 2-cycle for the relevant 
difference equation (for equations with a unique constant concentrated delay). 

5 Applications and Discussion 

The global attractivity results can be applied to any unimodal function / and any type of 
finite delay. For example, we can deduce the following result in the stream of [6], see also [8] 
and Proposition 2.1 in [25] . 

Theorem 9 Suppose (al)-(a5) hold, f is three times continuously differentiate and has the 
only critical point xq > (maximum) , 

^=jU-l{m) 2<0 f ° r x * x ° 

and \f'(K)\ < 1. Then any solution of (EP ; (G|) tends to K as t — > oo. 

As an application, consider the Nicholson's blowflies equation with a distributed delay 

x{t) = -Sx(t) +p f x{s)e- ax{s) d s R(t,s), (38) 
Jh(t) 

where (a2)-(a5) hold. 

If p < 5 then by Theorem [6] all positive solutions of (1581 go to extinction, i.e., x(t) — > 
p 

as t — > oo. If 1 < — < e 2 then by Theorem [5] all positive solutions of (BHD t en d to the 
o 

positive equilibrium K = — In (- ). However, since sustainable oscillations were observed in 

a \oJ 

experiments [10], we speculate that in the real problem the ratio p/5 was above e 2 and delays 
significant enough. 

Particular cases of ( |38l) are the equation with a variable delay 

x(t) = -5x{t) + px{h(t))e- ax{h{t)) , (39) 
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including the original equation (J4]) with a constant delay, the integrodifferential equation 

x(t) = -8x(t) + p [ k(t- s)x(s)e- ax{s) ds, (40) 
Jh(t) 



where 

IhU) 

and the mixed equation 



/ k{t - s) ds = 1 for any t > 0, (41) 
Jh(t) 



x(t) = -5x(t) + apx{g{t))e- axi9{t)) + (1 - a)p [ k(t - s)x(s)e~ ax{s) ds, < a < 1. (42) 

By Theorems El [6] and [7] we conclude the following result. 

Theorem 10 If p < 5 then all positive solutions of (3&) (and thus any of $MI\) , (TJty , §4]^) ) 

go to extinction, if 5 < p < 5e 2 , all positive solutions tend to the positive equilibrium. If 
p > e 2 and 

limsup(t -hit)) < ——, (43) 

<->oo p + o 

then all positive solutions of f3W tend to the positive equilibrium. 



Proof. The results for p < 5e 2 are immediate corollaries of Theorems [5] and [6j Since the 

p _ 

absolute value of the derivative of f{x) = -xe ax does not exceed p/S (which is attained at 
x = 0), then (JHTJ implies (1131. □ 



Remark 4 We can apply Remark[J\ to improve condition fcffi) . In fact, L* defined in 
does not exceed p/ (5e 2 ) (the absolute value of the minimum of the derivative which is attained 
at x = 2/a) if p > 5e 2 and the stability condition becomes 

limsup(t-/i(t)) < I (44) 
t^oo p/e z + d 

We also notice that for unimodal functions such that their equilibrium (the fixed point x = 
f(x)) exceeds maximum point x max and such that the minimum of the derivative is less 
than -1 (otherwise, under some additional conditions which are satisfied for the Nicholson's 
blowflies and the Mackey-Glass equations, the relevant difference equations are stable fEE/), 
then L* defined in p^ ) does not exceed the absolute value of the minimum of the derivative. 

fix) - K 

— — satisfies g(0) = 1, < g(x) < — min f'(x) in 

X — K xG[x max ,oo) 

[^max, oo) . The value in the right hand side of the latter inequality exceeds 1 and g is de- 
creasing in [0, x*\, where x* < x max is such point that f{x*) = K, while in [x*,x max ) function 
g{x) is less than g(x ma , x ), which does not exceed the absolute value of the minimum of the 
derivative. 



Really, function g(x) 
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Let us compare Theorem [10] and Remark [4] to some known global stability results for 
the Nicholson's blowflies equation. For example, [22J contains the following global stability 
condition 

-1) (f-l)<l (45) 

for equation (TJJ with a constant delay. 

Condition (jl5|) is slightly sharper for (T4]) than (j4"51) but fH3l) is applicable to equations 
with variable delays and integrodifferential equations as well. Condition (jUJ) is sharper than 
15]), see Fig. HI It is also applicable to equations with a variable and/or distributed delay. 
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Figure 4: The values of delay for which (jlj) is asymptotically stable are given for 5=1, 
p > 5e 2 > 7 and conditions f|44l) . (j4"3]) and fj45l) . respectively. 



For completeness of references, we remark that oscillation of fl38|) was studied in |26j . 
All the main results are also relevant for the Mackey-Glass equation with a distributed 
delay 

f ax(s) 

x(t) = I - ^ J d 8 R(t, s) - bx(t), a > 0, b > 0, 7 > 0, (46) 



h{t) 



1 + x"f(s) 



which also involves the equation 

aax(g(t)) 



x(t) = -bx(t) + — v^v^ + (1 - a)a [ kit - s) — r/.s. () < n <^ I . ( 17) 

where fT4T!) is satisfied, as well as its special cases of the equation with a variable delay 
(a = 1) and the integrodifferential equation (a = 0). If a < b then all solutions of (jlBTl go to 
extinction. If a > 6 and 

2a 



x(s) 



7 < 



a — b 



(4* 



then all solutions tend to the positive equilibrium = 
Again, Theorems [5J [6] and [7| imply the following result. 



1/7 
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Theorem 11 If a < b then all positive solutions of fijh] ) (and thus of ) go to extinction. 

76 

If either < 7 < 2 or 7 > 2 andb < a < , then all positive solutions tend to the positive 

7-2 

76 

equilibrium K. If ^ > 2, a > and 

7-2 

lim sup (t-h(t))< ,J - , , (49) 



then all positive solutions of (46) tend to the positive equilibrium. 

ax 

Proof. We remark that for 7 < 1 the function fix) = — is monotone increasing, 

K ' 6(1 + 

and we have stability for any delay. For 1 < 7 < 2 the relevant difference equation x n+ \ = 
f(x n ) is locally asymptotically stable, we apply Theorem It is known that for the relevant 
difference equation local asymptotic stability implies global asymptotic stability (see |25j, 
Example 3.2), thus the difference equation is globally asymptotically stable whenever 



\f'(K)\ 



b 

1 - 7 + -7 
a 



< 1, where 7 > 2, 



or a < 76/(7 — 2). According to Remark HI we can take the absolute value of the minimum 

7 + 1XV7 



of the derivative which is attained at ( ) as L*\ 

7 _1 



£ . _ q(7 - 1) 2 
467 

Application of Theorem [7] and Remark H] completes the proof. □ 
Finally, let us formulate some open problems. 

1. Generally, absolute stability results are incorrect for systems with infinite memory. 
Deduce stability results when (a2) is not satisfied but measure d s R(t, s) decays expo- 
nentially with memory 

ft-T 



/t — T 
d s R(t, s) < Se 
-00 



2. For arbitrary f(x) satisfying (al) and the equation with a constant delay, characterize 
attractive sets. 

3. Obtain stability results for equations with a distributed delay in the multistability 
case: there are several points satisfying f(x) = x. If Ki is the first positive fixed point 
and f(x)<x for < x < Ki, describe initial conditions leading to extinction (Alley 
effect). 
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4. Would the results of the present paper remain valid for the equation 



x(t) = r(t) 



f(f x(s)d s R(t,s)) -x(t) 
. \Jh(t) J 



,t>0, 



if (al)-(a5) hold? 
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